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ABSTRACT 

We investigate the asymmetry of the scattering cross section of radiation 
around Lya by atomic hydrogen, which may be apphed to analyses of scat- 
tering media with high column neutral hydrogen densities including damped 
Lya absorption systems of quasars. The exact scattering cross section is given 
by the Kramers-Heisenberg formula obtained from the fully quantum mechan- 
ical second-order time dependent theory, where, in the case of hydrogen, each 
matrix element is given in a closed analytical form. The asymmetric devia- 
tion of the scattering cross section from the Lorentzian near the line center is 
computed by expanding the Kramers-Heisenberg formula in terms of Au/uLya, 
where uiya is the angular frequency of the Lya transition and Auj is the devia- 
tion of incident radiation from ujiya- To the first order of Auj/ujiya, we obtain 
cr(cj) = crT(0.5/i2u;Lj^Q,/Au;)^(l — l.lS^Auj /ujLya), where ax is the Thomson scat- 
tering cross section and /12 = 0.4162 is the oscillator strength for the Lya transi- 
tion. With this deviation, the line center of the damped wing profile apparently 
shifts blueward of the true Lya line center. In the case of a damped Lya system 
with a H I column density 5 x 10^^ cm^^, the apparent line center shift relative 
to the true center amounts to 0.2 A resulting in an underestimation of redshift 
by Az ~ 10~^. A measurable underestimation by an amount of Az ~ 10^^ is 
expected for absorbing systems with Nhi > 4 x 10^^ cm~^. 

Subject headings: atomic physics — radiative transfer — cosmology — quasar 



1. Introduction 



The scattering process of electromagnetic radiation by an atom is regarded as a second- 
order quantum mechanical process involving the annihilation of an incident photon and the 
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creation of a scattered photon. Therefore, from the second-order time dependent pertur- 
bation theory, we can obtain the scattering cross section, which is known as the Kramers- 
Heisenberg formula (e.g. Sakurai 1969). The exact scattering cross section can be obtained 
after summing all bound state contributions and integrating over all continuum state contri- 
butions related with the dipole processes of annihilation and creation of photons. Hydrogen 
being a single electron atom, the energy eigenstate wavefunctions of hydrogen are analyti- 
cally known and hence each matrix element appearing in the Kramers- Heisenberg formula 
is given in a closed form. 

Near resonance, the scattering cross section is excellently approximated by the Lorcntzian 
characterized by the narrow radiation damping term that can be interpreted as the life time 
of the excited state. The Lorentzian dependence is obtained since the cross section is domi- 
nated by the single contribution from the excited state relevant to the resonance whereas the 
contributions from the other states are almost completely negligible. This Lorentzian ap- 
proximation is quite generally used in the astronomical community for analyses of resonance 
lines. 

Far from resonance, the radiation damping term is negligibly important and the scat- 
tering cross section is approximately proportional to Ao;"^, where Ao; is the difference of the 
angular frequency of the incident radiation and that of the Lya line center, ULyo,. This cross 
section given by the Lorentzian is symmetric with respect to the line center in frequency 
space. However, the true scattering cross section deviates from Acj^^ very far from the line 
center as the contributions increase from the other energy cigcnstatcs including both bound 
and continuum states. The zcroth order term cr oc {LUiya/ being symmetric with respect 
to the line center in frequency space, the first order term Au/uiLya for this deviation will 
contribute to the asymmetry of the scattering cross section. Therefore, a very careful quan- 
titative consideration of this effect should be taken in order to obtain an accurate estimation 
of the true line center of high H I column density absorption systems. 

H I absorption systems are easily found in typical quasar spectra and are very useful to 
investigate the baryonic matter content and the large scale structure in the universe. The 
quasar H I absorption systems include the Lya forest, Lyman limit systems and damped 
Lya systems, which are classified according to the H I column density N^^j. A number 
of damped Lya systems (DLA's) are known to exhibit H I column density in excess of 
10^^ cm~^ (e.g. Turnshek & Rao, 1998). Very high H I column density media may be found 
in the early universe, when the reionization is not complete and the universe is only partially 
ionized by the first objects (Gunn & Peterson 1965, Scheuer 1965). In these high column 
H I density media, the correct redshift of the absorbing system can be obtained after proper 
considerations of the asymmetric scattering cross section in frequency space. 
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In this paper, we calculate the asymmetric deviation of the scattering cross section of 
radiation around Lya by atomic hydrogen and apply our results to high H I column density 
scattering media that may be important in observational cosmology. 



2. Calculation 

2.1. The Kramers-Heisenberg Formula 

We start with the Kramers-Heisenberg formula that is obtained from the fully quantum 
mechanical second-order time dependent theory (e.g. Sakurai 1969, Mcrzbacher 1970). In 
terms of the matrix elements of the dipolc operator (ex), the differential cross section ^ 
for incident radiation with angular frequency uj is given by the Kramers-Heisenberg formula, 
which can be written as 
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where me is the electron mass, Tq = e^/mgC^ is the classical electron radius, (x • e(°'^)i„ 
represents the matrix element of the position operator between the ground Is state and 
excited np state, and e",e"' are polarization vectors of incident and scattered radiation, 
respectively. Here, uJni is the angular frequency between the bound np state and the ground 
Is state, and uOn'i is the angular frequency of the transition between Is and the continuum 
n'p state. In the atomic units adopted in this work, the bound np state has an energy 
eigenvalue En = and correspondingly ujni = (|)(1 — 1/n^). Similarly, En' = 2(^7j7 and 
a;„a = (|)(l + lK). 

In the limit uj oo, with the use of the closure relation it can be shown that the 
Kramers-Heisenberg formula yields a constant scattering cross section independent of the 
wavelength of incident radiation, which is the Thomson scattering cross section da/dQ — 
rgle^") • e^^'^p. In the opposite limit where u; <S 0021 — ojiya, another use of the completeness 
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of the intermediate states gives 



da 



+ (x-e("))u(x-e('^'On] 



(2) 



which shows the well-known lj^ dependence of the Rayleigh scattering cross section in the 
low energy limit. 

In order to investigate the behavior of the cross section around Lya, we introduce 

Acu = 00— L02\ (3) 

and expand the Kramers- Heisenberg formula in terms of Aiojio^x- 

In the first term inside the summation of Eq. 1 with n = 2, we can write 
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In the same term for n > 3, we may write 
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In a similar way, for n >2, the second term in the summation in Eq. 1 can be written as 
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Similar expressions for continuum states are obtained in a straightforward manner. 



Making use of the Wigner-Eckart theorem wc may separate the matrix elements of the 
position operator into the angular part and the radial part, where the radial part is given 
by the reduced matrix elements < f \\ x \\ i > . (e.g. Mcrzbacher 1970). The angular part 
gives the same dipolc type angular distribution as that for Thomson scattering or classical 
Rayleigh scattering (e.g. Ahn, Lee & Lee 2000, Lee 1997, Schmid 1989). 



Substituting Eqs 5,6 into Eq. 1, we obtain 
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where the atomic units have been used and angular integration has been performed. Here, the 
Thomson scattering cross section ar — Snr^/S — 0.665 x 10~^^ cm^. The matrix elements of 
the dipole operators < x >„i, < x >n',i are easily found in textbooks on quantum mechanics 
(e.g. Berestetski, Lifshitz & Pitaevskii 1971, Saslow & Mills 1969, Bethe & Salpeter 1967). 
The matrix elements are given by 



< X >ni—< ls\x\np >-- 
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for the bound states. Here, the Bohr radius = fi^/iTLe"^ is set equal to 1 in atomic units. 
For the continuum states, the corresponding values are given by 



< X >n'i=< ls\x\n'p >- 

where we also adopt the atomic units. 
Eq. 7 can be written as 
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where the coefficients Aj^ are obtained through the following relations 
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Here, /12 = 0.4162 is the oscillator strength for the Lya transition. 

The computation of the coefficients can be done in a straightforward way, from which 
we obtain a list of these coefficients up to A; = 5 

A^/Aq = -0.8961 

A2/AQ = -1.222 X 10^ 

= -5.252 X 10^ 

AJAq = -2.438 X 10^ 

A^/Aq = -1.210 X 10^ (12) 

Therefore, up to the first order in Aa;/a;2i, the scattering cross section of radiation 
around Lya 




In Fig. 1, we show the scattering cross section multiplied by a column density Nfji = 
4 X 10^^ cm~^, that is the scattering optical depth of a neutral medium with this column 
density. The solid line represents the result obtained from the Kramers-Heisenberg formula, 
the dotted line depicts the first order approximation to it obtained in Eq. 13 and the dashed 
line shows the zeroth approximation proportional to Acu"^ obtained from the Lorentzian by 
ignoring the damping part. The inclusion of the first order term implies the introduction of 
asymmetry of the scattering cross section with respect to the line center in frequency space. 
As is apparent in Fig. 1, the scattering cross section is smaller than the Lorentzian blueward 
of the line center, whereas it is larger in the opposite part. 

2.2. Applications to High Nhi Media 

In Fig. 2 we show the absorption profile e"'^^'^) that is expected from a damped Lya 
system with a column density Nhi = 5 x 10^^ cm~^. The thick sohd fine represents the 
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result obtained using the first-order corrected formula and the dotted fine is the absorption 
profile obtained using the Voigt profile by shifting the center wavelength by an amount 
AA = 0.187 A to the blue part. The wavelength shift by this amount is determined using the 
least square method that best fits the thick solid line. Because in the observational data it is 
difficult to obtain the exact continuum level, we minimize the difference |e^^^('^) — e"'^^''*^'^^ p in 
the wavelength interval between 1180 A and 1250 A, where the suppression of the continuum 
is larger than 20 percent. 

Therefore, if we apply the Voigt profile fitting process to observational data for an H 1 
scattering system with Nhi = 5 x 10^^ cm~^, we may underestimate the rcdshift of the H 1 
system by an amount of Az = 0.187/1215.67 = 1.54 x 10~^ or in terms of the velocity scale 
Av = 39 km s^^ 

Rao & Turnshek (2000) investigated the damped Lya absorption system along the line 
of sight to the quasar PKS 1127-145, which was best fitted with the Voigt profile with 
z — 0.3127 (see also Rao et al. 2003). Hence, the first order correction from our quantum 
mechanical consideration may affect the fourth decimal place of z, which is not critically 
important considering the accuracy of Az — 0.001 quoted by the authors. Among DLA's so 
far investigated, the highest H I column density is ~ N^hi — 5x 10^^ cm~^ (Rao & Turnshek 
2000). 

According to Braun (1997), gas rich galaxies may contain a medium with Nhi = 4 x 
10^^ cm~^. In Fig. 3 we show the scattering optical depth for this H I column density. For a 
higher column density system with Nhi = 4 x 10^^ cm~^, we obtain A A = 1.17 A, which may 
result in underestimation of the redshift by Az = 0.962 x 10~^. This amounts to a velocity 
shift of Av = 289 km s~^, and may affect significantly the exact redshift determination of 
the absorbing system. 

3. Discussion 

Gavrila (1967) provided a quantum mechanical computation of the Rayleigh scattering 
cross section using the Green function in momentum space. Ferland (2001) obtained the 
power fit to the result of Gavrila (1967) and incorporated it in his photoionization code 
'Cloudy,' in which the formula is useful for radiation with A > 1410 A and for radiation 
nearer the Lya line center Voigt function is proposed to be used. In the low energy limit 
where cu -C utiya, the scattering cross section tends to be proportional to u;^, which is also 
expected from classical physics. This approximation is quite useful in dealing with basic 
radiative transfer, because the Voigt function (or equivalently the Lorentzian function) well 
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approximates the true scattering cross section. 

Therefore, in this work, we emphasize the asymmetric deviation of the true scattering 
cross section from the Lorentzian in frequency space, which may result in underestimation 
of the exact wavelength of the Lya line center. Careful comparisons of redshifts determined 
using other lines may provide important clues to the kinematics of the absorbing systems. In 
particular, Prochaska & Wolfe (1997, 1998) proposed that the absorption profiles observed in 
many DLA's are consistent with the hypothesis that they are formed in a rotating disk. The 
apparent line shift due to the asymmetric deviation around Lya may affect the kinematic 
analysis of rotating galaxies. 

Scattering of Lya photons is also important in search for the first objects that are 
responsible for the reionization of the universe (Gunn & Peterson 1965, Scheuer 1965). In 
several high red shift quasars with z > 6.2 the presence of the Gunn-Peterson trough has 
been reported blueward of Lya (Becker et al. 2001, Fan et al. 2003). The Gunn-Peterson 
optical depth will be important in the red wing part of Lya (e.g. Miralda-Escude 1998). In 
Eq. (13) the correction from the Lorentzian cross section is about 1 % when we consider a 
width of the damped profile of about 2000 km s~^, for the characteristic expected value of 
the intergalactic medium (IGM) density at around z — 6. Hence, the change in the cross 
section is equivalent to a change in the redshift of the Gunn-Peterson trough of about 1 % 
of th ewidth of the damped profile or 20 km s~^, when only the first order term in Au}/u}2i 
is considered. 

Another high column density media may be found in the atmospheres of supergiants, 
in which classical Rayleigh scattering is important (e.g. Isliker Nussbaumer & Vogel 1989). 
Also about a half of the symbiotic stars believed to consist of a heavily mass losing giant 
and a white dwarf exhibit Raman scattered O VI features around 6830 A and 7088 A 
(Schmid 1989, Nussbaumer, Schmid, & Vogel 1989). These features are originated from the 
resonance doublet O VI 1032, 1038 line photons that are scattered in a neutral region with 
Nhi ~ 10^^ cm~^. However, due to resonant scattering by neutral hydrogen in interstellar 
space, it will be difficult to accurately determine the asymmetry of the scattering cross 
section. 

The author is very grateful to an anonymous referee for useful comments. This work 
is a result of research activities of the Astrophysical Research Center for the Structure 
and Evolution of the Cosmos (ARCSEC) funded by the Korea Science and Engineering 
Foundation. 
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Fig. 1. — The scattering cross section around Lya. The result obtained from the Kramers- 
Heisenberg formula is shown by the solid line, the first order approximation to it obtained 
in Eq. 13 and the zeroth approximation that is proportional to Ao;"^ obtained from the 
Lorentzian by ignoring the damping part. The scattering cross section is smaller than the 
Lorentzian blueward of the line center, whereas it is larger in the opposite part. 
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Fig. 2. — The absorption profile e"^*-'^-' expected from a damped Ly« system witli a column 
density = 5 x 10^^ cm~^. The thick solid line represents the result obtained using the 
first-order corrected formula and the dotted line is the absorption profile obtained using the 
Voigt profile by shifting the center wavelength by an amount AA = 0.187 A to the blue part. 
The wavelength shift by this amount is determined using the least square method that best 
fits the thick solid line. 
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Fig. 3. — The absorption profile e^^*^'*') expected from a damped Lya system witli a column 
density Nhi = 4 x 10^^ cni^^. The thick solid line represents the first-order correction. 
The best-fit wavelength shift from the Lorentzian fit is A A = 1.17 A, which may result in 
underestimation of the redshift by Az — 0.962 x 10~^. 



